Event-by-event simulation of quantum phenomena* 
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In this talk, I discuss recent progress in the development of simulation algorithms that do not 
rely on any concept of quantum theory but are nevertheless capable of reproducing the averages 
computed from quantum theory through an event-by-event simulation. The simulation approach is 
illustrated by applications to single-photon Mach-Zehnder interferometer experiments and Einstein- 
Podolsky-Rosen-Bohm experiments with photons. 
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I. INTRODUCTION 

Computer simulation is widely regarded as comple- 
mentary to theory and experiment The standard 
procedure is to start from one or more basic equations of 
physics and to apply existing or invent new algorithms to 
solve these equations. This approach has been highly suc- 
cessful for a wide variety of problems in science and en- 
gineering. However, there are a number of physics prob- 
lems, very fundamental ones, for which this approach 
fails, simply because there are no basic equations to start 
from. 

Indeed, as is well-known from the early days in the de- 
velopment of quantum theory, quantum theory has noth- 
ing to say about individual events @, H, 0| ■ Reconciling 
the mathematical formalism that does not describe indi- 
vidual events with the experimental fact that each ob- 
servation yields a definite outcome is referred to as the 
quantum measurement paradox and is the most funda- 
mental problem in the foundation of quantum theory Q . 

If computer simulation is indeed a third methodology, 
it should be possible to simulate quantum phenomena on 
an event-by-event basis. For instance, it should be possi- 
ble to simulate that we can see, with our own eyes, how 
in a two-slit experiment with single electrons, an inter- 
ference pattern appears after a considerable number of 
individual events have been recorded by the detector || . 

In view of the quantum measurement paradox, it is 
unlikely that we can find such a simulation algorithm 
by limiting our thinking to the framework of quantum 
theory. Of course, we could simply use pseudo-random 
numbers to generate events according to the probabil- 
ity distribution that is obtained by solving the time- 
independent Schrodinger equation. However, that is not 
what we mean when we say that within the framework of 
quantum theory, there is little hope to find an algorithm 
that simulates the individual events and reproduces the 
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expectation values obtained from quantum theory. The 
challenge is to find algorithms that simulate, event-by- 
event, the experimental observations that, for instance, 
interference patterns appear only after a considerable 
number of individual events have been recorded by the 
detector without first solving the Schrodinger equa- 

tion. 



In a number of recent papers 0, i, S [M EH, E E 0, 
[l5j . we have demonstrated that locally-connected net- 
works of processing units can simulate cvcnt-by-cvcnt, 
the single-photon beam splitter and Mach-Zchnder in- 
terferometer experiments of Grangier et al. [□]. Further- 
more, we have shown that this approach can be gener- 
alized to simulate universal quantum computation by an 
event- by-event process [1, [l(| [H[ , and that it can be used 
to simulate real Einstein-Podolsky-Rosen-Bohm (EPRB) 
experiments 0, [3 EBl • Therefore, at least in principle, 
our approach can be used to simulate all wave interfer- 
ence phenomena and many-body quantum systems using 
particle-like processes only. Our work suggests that we 
may have discovered a procedure to simulate quantum 
phenomena using event-based processes that satisfy Ein- 
stein's criterion of local causality. 



This talk is not about interpretations or extensions of 
quantum theory. The fact that there exist simulation al- 
gorithms that reproduce the results of quantum theory 
has no direct implications on the foundations of quantum 
theory: These algorithms describe the process of generat- 
ing events on a level of detail about which quantum the- 
ory has nothing to say [!, 3] . The average properties of 
the data may be in perfect agreement with quantum the- 
ory but the algorithms that generate such data are out- 
side of the scope of what quantum theory can describe. 
This may sound a little strange but it may not be that 
strange if one recognizes that probability theory does not 
contain nor provides an algorithm to generate the values 
of the random variables either, which in a sense, is at the 
heart of the quantum measurement paradox [lfj . 
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Single-photon Mach-Zehnder interferometer 
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FIG. 1: (color online) Snapshot of an interactive event-by-event simulator of a Mach-Zehnder interferometer jig ]. The main 
panel shows the layout of the interferometer. Particles emerge from a source (not shown) located at the bottom of the left-most 
vertical line. After leaving the first beam splitter in either the vertical or horizontal direction, the particles experience time 
delays that are specified by the controls on the lines. In this example, the time delays correspond to the phase shifts cj>o = 35° 
and (pi — 322° in the wave mechanical description. The thin, 45°-tilted lines act as perfect mirrors. When a particle leaves 
the system at the top right, it adds to the count of either detector N2 or A3. Additional detectors (No, iVi) count the number 
of particles on the corresponding lines. The other cells give the ratio of the detector counts to the total number of particles 
(messages) processed and also the corresponding probability of the quantum mechanical description. At any time, the user can 
choose between a strictly deterministic and a stochastic event-by-event simulation by pressing the buttons at the top of the 
control panel. 



II. SINGLE-PHOTON MACH-ZEHNDER 
INTERFEROMETER 

Figure Q] shows the schematic diagram of a Mach- 
Zehnder interferometer [ijj ■ From Maxwell's theory of 
classical electrodynamics it follows that the intensity of 
light recorded by detectors N2 and N3 is proportional to 
cos 2 <p/2 and sin 2 <p/2, respectively [171] ■ Here 4> = 4>i — 02 
is the phase difference that expresses the fact that de- 
pending on which path the light takes to travel from the 
first to the second beam splitter, the optical path length 
may be different [l7| . 

It is an experimental fact that when the Mach-Zehnder 
interferometer experiment is carried out with one photon 
at a time, the number of individual photons recorded 
by detectors N2 and N3 is proportional to cos 2 0/2 and 
sin 2 (f>/2 Q, in agreement with classical electrodynam- 



ics. In quantum physics [l8j . single-photon experiments 
with one beam splitter provide direct evidence for the 
particle-like behavior of photons. The wave mechanical 
character appears when one performs interference exper- 
iments with individual particles 3, 6]. Quantum physics 
"solves" this logical contradiction by introducing the con- 
cept of particle- wave duality [3J. 

In this section, we describe a system that does not 
build on any concept of quantum theory yet displays the 
same interference patterns as those observed in single- 
photon Mach-Zehnder interferometer experiments [6(. 
The basic idea is to describe (quantum) processes in 
terms of events, messages, and units that process these 
events and messages. In the experiments of Grangier et 
al. [f| , the photon carries the message (a phase) , an event 
is the arrival of a photon at one of the input ports of a 
beam splitter, and the beam splitters are the process- 
ing units. In experiments with single photons, there is 



no way other than through magic, by which a photon 
can communicate directly with another photon. Thus, 
it is not difficult to imagine that if we want a system 
to exhibit some kind of interference, the communication 
among successive photons should take place in the beam 
splitters. 

In this talk, we consider the simplest processing unit 
that is adequate for our purpose, namely a standard lin- 
ear adaptive filter The processing unit receives a 
message through one of its input ports, processes the 
message according to some rule (see later), and sends 
a message (carried by the messenger, that is a photon) 
through an output port that it selects using a pseudo- 
random number, drawn from a distribution that is deter- 
mined by the current state of the processing unit. Other, 
more complicated processing units that operate in a fully 
deterministic manner are described elsewhere 
though the sequence of events that the different types of 
processing units produce can be very different, the quan- 
tities that are described by quantum theory, namely the 
averages, are the same. The essential feature of all these 
processing units is their ability to learn from the events 
they process. Processing units that operate according to 
this principle will be referred to as deterministic learning 
machines (DLMs) @,G3. 

By connecting an output channel to the input channel 
of another DLM, we can build networks of DLMs. As the 
input of a network receives an event, the corresponding 
message is routed through the network while it is being 
processed and eventually a message appears at one of the 
outputs. At any given time during the processing, there 
is only one output-to-input connection in the network 
that is actually carrying a message. The DLMs process 
the messages in a sequential manner and communicate 
with each other by message passing. There is no other 
form of communication between different DLMs. The 
parts of the processing units and network map one-to- 
one on the physical parts of the experimental setup and 
only simple geometry is used to construct the simulation 
algorithm. Although networks of DLMs can be viewed as 
networks that are capable of unsupervised learning, they 
have little in common with neural networks. It obvious 
that this simulation approach satisfies Einstein's criteria 
of realism and local causality 0] . 



A. Beam splitter 

Figure [2] shows the schematic diagram of a DLM 
that simulates a beam splitter, event-by-event. We la- 
bel events by a subscript n > 0. At the (n + l)th 
event, the DLM receives a message on either input chan- 
nel or 1, never on both channel simultaneously. Ev- 
ery message consists of a two-dimensional unit vector 
y n+ i — (yi,n+i, 2/2, n+i)- This vector represents the phase 
of the event that occurs on channel (1). Although it 
would be sufficient to use the phase itself as the message, 
in practice it is more convenient to work with the cosine 
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FIG. 2: Diagram of a DLM that performs an event-by-event 
simulation of a single-photon beam splitter (BS) The 
solid lines represent the input and output channels of the BS. 
Dashed lines indicate the flow of data within the BS. 



(yi.n+i) and sine (j/2,n+i) of the phase. 

The first stage of the DLM (see Fig. [5]) stores the mes- 
sage y n +i in its internal register Y&. Here, k = (1) if 
the event occurred on channel (1). The first stage also 
has an internal two-dimensional vector x = (x , Xi) with 
the additional constraints that Xi > for i = 0, 1 and 
that xq + X\ = 1. As we only have two input channels, 
the latter constraint can be used to recover X\ from the 
value of xo ■ We prefer to work with internal vectors that 
have as many elements as there are input channels. After 
receiving the (n + l)-th event on input channel k = 0, 1 
the internal vector is updated according to the rule 

%i,n+i = axi^n + 1 - a if i = k, 

Xi, n +i = ax it n if i ^ k, (1) 

where < a < 1 is a parameter. By construction 
Xi.n+i > for i = 0, 1 and xo. n +i + a^i.n+i = 1- Hence 
the update rule Eqs. ([T]) preserves the constraints on the 
internal vector. Obviously, these constraints are neces- 
sary if we want to interpret the Xk, n as (an estimate of) 
the probability for the occurrence of an event of type k. 

The second stage of the DLM takes as input the values 
stored in the registers Yq, Yj., x and transforms this data 
according to the rule 
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where we have omitted the event label (n+1) to simplify 
the notation. Note that the second subscript of the Y- 
register refers to the type of input event. 

The third stage of the DLM in Fig. [2] responds to the 
input event by sending a message w n+ i = (Yo^oy/xo — 
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Yi,i ^0,1 V^i + Yifiy/%o)/ V2 through output chan- 
nel if Wq n+1 + w\ n+1 > r where < r < 1 is a uniform 
random number. Otherwise the back-end sends the mes- 
sage z„ + i = (Fo^V^ _ ^bOv / ^'^o,ov / ^o + ^i,iv / ^i)/v^ 
through output channel 1 . Finally, for reasons of internal 
consistency of the simulation method, it is necessary to 
replace w n+ i by w„ + i/||w„ + i|| or z n+1 by z„ + i/||z n+ i || 
such that the output message is represented by a unit 
vector. 

It is almost trivial to perform a computer simulation of 
the DLM model of the beam splitter and convince oneself 
that it reproduces all the results of quantum theory for 
this device @. With only a little more effort, it can be 
shown that the input-output behavior of the DLM is, on 
average, the same as that of the (ideal) beam splitter. 

According to quantum theory, the probability ampli- 
tudes (bo,bi) of the photons in the output modes and 
1 of a beam splitter (see Fig. ^ are given by @, [l^, [2(| 

( b ° \ = _L ( a ° + iCl1 \ = _L ( 1 i \ ( a o\ fo\ 

\bi J ^{a.+iaoj ^ \i 1 J \ *i J ' W 

where the presence of photons in the input modes or 1 is 
rep resented by the probability amplitudes (arj, ai) [rUlH 
|2G | . From Eq. [3l it follows that the intensities recorded 
by detectors No and Ni is given by 

,2 _ 1 + 2y/p (l - po) sm(t/j - V>i) 



| fci |2 _ 1 - 2 Vpo( 1 - Po) sin(-0 o - ^ 

On the other hand, the formal solution of Eq. JT]) reads 

n-i 

x„ = a"x + (l-a)^a"~ 1 -V J+ i, (5) 

i=0 

where x„ = (io,n,ii.n), and xo denotes the initial value 
of the internal vector. The input events are represented 
by the vectors v n+ i = (1,0) T or v„ + i = (0, 1) T if the 
n + 1-th event occurred on channel or 1 , respectively. 
Let po ((1 — po)) be the probability for an input event 
of type (1). Taking the average of Eq.© over many 
events and using < a < 1, we find that for large n, 
Xn ~ (po, 1 — Po) T ■ Therefore the first stage of the DLM 
"learns" the probabilities for events and 1 by processing 
these events in a sequential manner. The parameter < 
a < 1 controls the learning process. 

Using two complex numbers instead of four real 
numbers that enter Eq. @, identification of with 

^o.ov^o + iii,ov^o and a x with F ,iV^i~ + ^l.iV^T 
shows that the transformation stage plays the role of the 
matrix- vector multiplication in Eq.Q. By construction, 
the output stage receives as input the four real num- 
bers that correspond to b and b\. Thus, after the DLM 
has reached the stationary state, it will distribute events 
over its output channels according to Eq.(d]). Of course, 
this reasoning is firmly supported by extensive simula- 
tions 




§ [degrees] 

FIG. 3: Simulation results for the DLM-network shown in 
Fig- [U Input channel receives (cos^o, sin^o) with proba- 
bility one. A uniform random number in the range [0, 360] 
is used to choose the angle ipo- Input channel 1 receives no 
events. Each data point represents 10000 events (No + Ni — 
N 2 + N s = 10000). Initially the rotation angle <j> = and 
after each set of 10000 events, <j>o is increased by 10° . Markers 
give the simulation results for the normalized intensities as a 
function of <f> = <f>o — <f>\. Open squares: Nq/{Nq + Ni); Solid 
squares: ^2/(^2+^3) for <pi = 0; Open circles: ^2/(^2+^3) 
for 0i = 30°; Bullets: ^2/(^2 + ^3) for 4>i = 240°; Asterisks: 
^3/(^2 + N 3 ) for 0i = 0; Solid triangles: 7V 3 /(jV 2 + JV 3 ) for 
0i = 300°. Lines represent the results of quantum theory [18|. 



One may wonder what learning machines have to do 
with the (wave) mechanical models that we are accus- 
tomed to in physics. First, one should keep in mind 
that the approach that I describe in this talk is capa- 
ble of giving a rational, logically consistent description of 
event-based phenomena that cannot be incorporated in a 
wave mechanical theory without adding logically incom- 
patible concepts such as the wave function collapse [3[. 
Second, the fact that a mechanical system has some kind 
of memory is not strange at all. For instance, a pulse of 
light that impinges on a beam splitter induces a polar- 
ization in the active part (usually a thin layer of metal) 
of the beam splitter [13]. Assuming a linear response 
(as is usually done in classical electrodynamics), we have 
P(r,i) = x(r,i) * E(r, t) where "*" is a shorthand for 
the convolution. If the susceptibility x( r ;^) has a non- 
trivial time dependence (as in the Lorentz model [l7| for 
instance) , the polarization will exhibit "memory" effects 
and will "learn" from subsequent pulses. DLMs mimic 
this behavior in the most simple manner (see the convo- 
lution in Eq. ([5])). on an event-by-event basis. 

B. Mach-Zehnder interferometer 

Using the DLM of Fig. [5] as a module that simulates a 
beam splitter, we build the Mach-Zehnder interferometer 
by connecting two DLMs, as shown in Fig. [T] The length 
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FIG. 4: (color online) Schematic diagram of an EPRB experiment with photons [2lj . 



of each path from the first to the second beam splitter is 
made variable, as indicated by the controls on the hor- 
izontal lines. The thin, 45°-tilted lines act as perfect 
mirrors. 

In quantum theory, the presence of photons in the in- 
put modes or 1 of the interferometer is represented by 
the probability amplitudes (oto,ai) [HI]. The amplitudes 
to observe a photon in the output modes and 1 of the 
Mach-Zehnder interferometer (see Fig. [I} are given by 

(£)-(! D(? ,6) 

where bo and b\ are given by Eq. ([3]). In Eq. the 
entries e % ^ i for j = 0, 1 implement the phase shifts that 
result from the time delays on the corresponding path 
(including the phase shifts due to the presence of the 
perfect mirrors). 

C. Simulation results 

The snapshot in Fig. Q] is taken after N — 3030 parti- 
cles have been generated by the source. The numbers in 
the various corresponding fields clearly show that even af- 
ter a modest number of events, this event-by-event simu- 
lation reproduces the quantum mechanical probabilities. 
Of course, this single snapshot is not a proof that the 
event-by-event simulation also works for other choices of 



the time delays. More extensive simulations, an exam- 
ple of a set of results being shown in Fig. [3l demonstrate 
that DLM-networks accurately reproduce the probabil- 
ities ofquantum theory for these single-photon experi- 
ments 0, in, Gl El Eg. 



III. EPRB EXPERIMENTS 

In Fig. 0] we show a schematic diagram of an EPRB 
experiment with photons (see also Fig. 2 in [2l|). The 
source emits pairs of photons. Each photon of a pair 
propagates to an observation station in which it is ma- 
nipulated and detected. The two stations are separated 
spatially and temporally [2l| . This arrangement prevents 
the observation at station 1 (2) to have a causal effect on 
the data registered at station 2(f) [2l|. As the photon ar- 
rives at station i = 1, 2, it passes through an electro-optic 
modulator that rotates the polarization of the photon by 
an angle depending on the voltage applied to the modu- 
lator. These voltages are controlled by two independent 
binary random number generators. As the photon leaves 
the polarizer, it generates a signal in one of the two de- 
tectors. The station's clock assigns a time-tag to each 
generated signal. Effectively, this procedure discretizes 
time in intervals of a width that is determined by the 
time-tag resolution r [2lj . In the experiment, the firing 
of a detector is regarded as an event. 
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As we wish to demonstrate that it is possible to re- 
produce the results of quantum theory (which implicitly 
assumes idealized conditions) for the EPRB gedanken 
experiment by an event-based simulation algorithm, it 
would be logically inconsistent to "recover" the results of 
the former by simulating nonideal experiments. There- 
fore, we consider ideal experiments only, meaning that 
we assume that detectors operate with 100% efficiency, 
clocks remain synchronized forever, the "fair sampling" 
assumption is satisfied [22j, and so on. We assume that 
the two stations are separated spatially and temporally 
such that the manipulation and observation at station 1 
(2) cannot have a causal effect on the data registered 
at station 2 (1). Furthermore, to realize the EPRB 
gedanken experiment on the computer, we assume that 
the orientation of each electro-optic modulator can be 
changed at will, at any time. Although these conditions 
are very difficult to satisfy in real experiments, they are 
trivially realized in computer experiments. 

In general, on logical grounds (without counterfactual 
reasoning), it is impossible to make a statement about 
the directions of the polarization of particles emitted by 
the source unless we have performed an experiment to de- 
termine these directions. However, in a computer exper- 
iment we have perfect control and we can select any di- 
rection that we like. Conceptually, there are two extreme 
cases. In the first case, we assume that we know nothing 
about the direction of the polarization. We mimic this 
situation by using pseudo-random numbers to select the 
initial polarization. This is the case that is typical for a 
real EPRB experiment. In the second case, we assume 
that we know that the polarizations are fixed (but are 
not necessarily the same) , mimicking a source that emits 
polarized photons. A simulation algorithm that aims to 
reproduce all the results of quantum theory should be 
able to reproduce all these results for both cases without 
any change to the simulation algorithm except for the 
part that simulates the source [la , [Til [l5| . 

In the experiment, the firing of a detector is regarded 
as an event. At the nth event, the data recorded on 
a hard disk at station i = 1,2 consists of x n ^ = ±1, 



specifying which of the two detectors fired, the time tag 
t nt i indicating the time at which a detector fired, and 
the two-dimensional unit vector a n ^ that represents the 
rotation of the polarization by the electro-optic polarizer. 
Hence, the set of data collected at station i — 1, 2 during 
a run of N events may be written as 

T 4 = {x n<i = ±l,t nt i,a n< i\n = 1,...,N}. (7) 

In the (computer) experiment, the data {Ti,T2} may 
be analyzed long after the data has been collected [2l| . 
Coincidences are identified by comparing the time dif- 
ferences — tn,2\n = 1 7 ■ • • , N} with a time window 
W [2l|. Introducing the symbol to indicate that the 
sum has to be taken over all events that satisfy = a„.; 
for i — 1,2, for each pair of directions ai and a 2 of 
the electro-optic modulators, the number of coincidences 
C xy = C xy (ai,a 2 ) between detectors D Xi i (x = ±1) at 
station 1 and detectors D y2 {v — ±1) at station 2 is 
given by 

N 

C xy = Sx ^Jv,x n , a Q(W - \t nA -t n , 2 \), (8) 

n=l 

where O(i) is the Heaviside step function. We emphasize 
that we count all events that, according to the same cri- 
terion as the one employed in experiment, correspond to 
the detection of pairs. The average single-particle counts 
are defined by 

^ , X El,J=±l XC X y 

and 

E 2 (a 1 ,a 2 ) = — , (9) 

2-^ x , y =±l ^xy 

where the denominator is the sum of all coincidences. 

According to standard terminology, the correlation be- 
tween x — ±1 and y — ±1 events is defined by [23| 



Ex, H XyCxy xCyy Y.x,y v c j 



P(ai,a 2 ) 



Y c x „ v 

^-^>x,y ■*- y £-^ix 



(10) 



The correlation p(ai,a2) is +1 (—1) in the case that 
x = y (x = — y) with certainty. If the values of x and y 
are independent, the correlation p(ai,a2) is zero. Note 
that in general, the converse is not necessarily true but 
in the special case of dichotomic variables x and y, the 
converse is true [24j . 

In the case of dichotomic variables x and y, the cor- 



relation p(ai,a2) is entirely determined by the average 
single-particle counts Eq. © and the two-particle aver- 
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FIG. 5: (color online) Smax as a function of the time window 
W, computed from the data sets contained in the archives 
Alice.zip and Bob. zip that can be downloaded from Ref. [25]. 
Bullets (red): Data obtained by using the relative time shift 
A = 4 ns that maximizes the number of coincidences. Crosses 
(blue): Raw data (A = 0). Dashed line at 2\/2: S m ax if the 
system is described by quantum theory (see Section IIII C|l . 
Dashed line at 2: Smax if the system is described by the class 
of models introduced by Bell [2o |. 
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£7(ai,a 2 ) = 



Yl x ,y x VCxy 
^2x.v Cxy 



(11) 



C++ + C—+C+- + C-+ 

For later use, it is expedient to introduce the function 

5(a, b, c, d) = E{&, c) - E(a, d) + E(b, c) + E(b, d), 

(12) 

and its maximum 



max S(sl, b, c, d) 

a.b.c.d 



(13) 



In general, the values for the average single- 
particle counts Ei(a.i,&2) and £? 2 (ai,a 2 ) the coinci- 
dences C xy {sL\ 1 a 2 ), the two-particle averages £?(ai,a 2 ), 
S(a, b, c, d), and S max not only depend on the directions 
ai and a 2 but also on the time-tag resolution r and the 
time window W used to identify the coincidences. 



A. Analysis of real experimental data 

We illustrate the procedure of data analysis and the 
importance of the choice of the time window W by an- 
alyzing a data set (the archives Alice.zip and Bob. zip) 
of an EPRB experiment with photons that is publicly 
available (2f|. 

In the real experiment, the number of events detected 
at station 1 is unlikely to be the same as the number of 
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FIG. 6: (color online) Same as Fig. [5] except for the range of 
W. Bullets (red): Data obtained by using the relative time 
shift A = 4 ns that maximizes the number of coincidences. 
The maximum value of S m „ w 2.73 is found at W = 2 ns. 
Crosses (blue): Raw data A = 0. The maximum value of 
Smax ~ 2.89 is found at W = 3 ns. 



events detected at station 2. In fact, the data sets of 
Ref. show that station 1 (Alice.zip) recorded 388455 
events while station 2 (Bob. zip) recorded 302271 events. 
Furthermore, in the real EPRB experiment, there may be 
an unknown shift A (assumed to be constant during the 
experiment) between the times t n> i gathered at station 1 
and the times t n 2 recorded at station 2. Therefore, there 
is some extra ambiguity in matching the data of station 
1 to the data of station 2. 

A simple data processing procedure that resolves this 
ambiguity consists of two steps 27]. First, we make a 
histogram of the time differences t n> i — t m 2 with a small 
but reasonable resolution (we used 0.5 ns). Then, we fix 
the value of the time-shift A by searching for the time 
difference for which the histogram reaches its maximum, 
that is we maximize the number of coincidences by a 
suitable choice of A. For the case at hand, we find A = 
4 ns. Finally, we compute the coincidences, the two- 
particle average, and S max using the expressions given 
earlier. The average times between two detection events 
is 2.5 ms and 3.3 ms for Alice and Bob, respectively. The 
number of coincidences (with double counts removed) is 
13975 and 2899 for (A = 4 ns, W = 2 ns) and (A = , 
W = 3 ns) respectively. 

In Figs. [5] and [6] we present the results for S max as a 
function of the time window W. First, it is clear that 
Smax decreases significantly as W increases but it is also 
clear that as W — > 0, S ma x is not very sensitive to the 
choice of W Second, the procedure of maximizing 

the coincidence count by varying A reduces the maxi- 
mum value of S m ax from a value 2.89 that considerably 
exceeds the maximum for the quantum system (2y/2, see 
Section Till C| to a value 2.73 that violates the Bell in- 
equality (S max < 2, see Ref. [26h and is less than the 
maximum for the quantum system. 



8 



c 

=3 
O 

o 

<D 
o 
c 

CD 
"O 
CJ 

c 
'o 

CJ 

"O 

CD 
N 

IB 
E 



1 


1 1 rX tr 1 1 




0.8 


^ / I 




0.6 


* i \ 




0.4 






0.2 


* / l 
/ ^ \ 






x' J \ 




> 


<H«^X— X x 



-4-3-2-101234 
tp,l-tn,2 

FIG. 7: (color online) Normalized coincidence counts as a 
function of time tag difference t„ t i — t n< 2, computed from the 
data sets contained in the archives Alice.zip and Bob. zip [2^ |. 
using the relative time shift A = 4 ns that maximizes the 
number of coincidences. Bullets (red): 9i = and 62 = tt/8; 
Crosses (blue): 9i = and #2 = 37r/8. 



The fact that the "uncorrected" data (A = 0) violate 
the rigorous bound for the quantum system should not 
been taken as evidence that quantum theory is "wrong" : 
It merely indicates that the way in which the data of 
the two stations has been grouped in two-particle events 
is not optimal. There is no reason why a correlation 
between similar but otherwise unrelated data should be 
described by quantum theory. 

Finally, we use the experimental data to show that 
the time delays depend on the orientation of the polar- 
izer. To this end, we select all coincidences between -D+,i 
and -D +i 2 (see Fig. 2]) and make a histogram of the coin- 
cidence counts as a function of the time-tag difference, 
for fixed orientation 9\ — and the two orientations 
82 = 7r/8, Stt/8 (other combinations give similar results). 
The results of this analysis are shown in Fig. [7] The max- 
imum of the distribution shifts by approximately 1 ns as 
the polarizer at station 2 is rotated by 7r/4, a demonstra- 
tion that the time-tag data is sensitive to the orientation 
of the polarizer at station 2. A similar distribution of 
time-delays (of about the same width) was also observed 
in a much older experimental realization of the EPRB 
experiment [28|. 

According to Maxwell's equation, the birefringent 
properties of the optically anisotropic materials that are 
used to fabricate the optical elements (polarizers and 
electro-optic modulators), cause plane waves with differ- 
ent polarization to propagate with different phase veloc- 
ity 13 ) suggesting a possible mechanism for the time 
delays observed in experiments. As light is supposed 
to consist of non-interacting photons, this suggests, but 
does not prove, that individual photons experience a time 
delay as they pass through the electro-optic modulators 
or polarizers. Of course, strictly speaking, we cannot de- 



rive the time delay from classical electrodynamics: The 
concept of a photon has no place in Maxwell's theory. 
A more detailed understanding of the time delay mech- 
anism first requires dedicated, single-photon retardation 
measurements for these specific optical elements. 



B. Role of the coincidence window W 

The crucial point is that in any real EPR-type exper- 
iment, it is necessary to have an operational procedure 
to decide if the two detection events correspond to the 
observation of one two-particle system or to the obser- 
vation of two single-particle systems. In standard "hid- 
den variable" treatments of the EPR gedanken experi- 
ment 26], the operational definition of "observation of 
a single two-particle system" is missing. In EPRB-type 
experiments, this decision is taken on the basis of coinci- 
dence in time [H [H ^ . 

Our analysis of the experimental data shows beyond 
doubt that a model which aims to describe real EPRB 
experiments should include the time window W and that 
the interesting regime is W — ► 0, not W — > 00 as is 
assumed in all textbook treatments of the EPRB exper- 
iment. Indeed, in quantum mechanics textbooks it is 
standard to assume that an EPRB experiment measures 
the correlation [26| 



N 

(7(00) X Sj 

n=l 



(14) 



which we obtain from Eq. ([5]) by taking the limit W — > 
00. Although this limit defines a valid theoretical model, 
there is no reason why this model should have any bear- 
ing on the real experiments, in particular because exper- 
iments pay considerable attention to the choice of W . A 
rational argument that might justify taking this limit is 
the hypothesis that for ideal experiments, the value of 
W should not matter. However, in experiments a lot of 
effort is made to reduce (not increase) W [2l], [27| • 

As we will see later, using our model it is relatively 
easy to reproduce the experimental facts and the results 
of quantum theory if we consider the limit W — ► 0. Fur- 
thermore, keeping W arbitrary does not render the math- 
ematics more complicated so there really is no point of 
studying the simplified model defined by Eq. (fT4"|) : We 
may always consider the limiting case W — > 00 after- 
wards. 



C. Quantum theory 

According to the axioms of quantum theory re- 
peated measurements on the two-spin system described 
by the density matrix p yield statistical estimates for the 
single-spin expectation values 



Si (a) 



(cri • a) 



E 2 (h) 



(<T2 ■ b) 



(15) 
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and the two-spin expectation value 

£(a,b) = (cri-acr 2 -b), (16) 

where Ci = (erf, of, of) are the Pauli spin- 1/2 matrices 
describing the spin of particle i = 1,2 jj], and a and b 
are unit vectors. We have introduced the tilde to dis- 
tinguish the quantum theoretical results from the results 
obtained from the data sets {Ti,T2}. The state of a 
quantum system of two S = 1/2 objects is completely de- 
termined if we know the expectation values -E'i(a), ^(b), 
and E(a, b). 

It can be shown that |5(a, b, c, d)| < 2V2 [3(|, in- 
dependent of the choice of p. If the density matrix 
p = Pi® P2 factorizes (here pi is the 2x2 density matrix 
of spin i), then it is easy to prove that |5(a, b, c, d)| < 
2. In other words, if max a b,c,d S(a, b, c, d) > 2, then 
P 7^ Pi ® Pi i and the quantum system is in an entangled 
state. Specializing to the case of the photon polarization, 
the unit vectors a, b, c, and d lie in the same plane and 
we may use the angles a, a' , 0, and (3' to specify their 
direction. 

The quantum theoretical description of the EPRB ex- 
periment assumes that the system is represented by the 
singlet state |*) = {\H) X \V) 2 - \V) X \H) 2 ) /V2 of two 
spin-1/2 particles, where H and V denote the horizon- 
tal and vertical polarization and the subscripts refer 
to photon 1 and 2, respectively. For the singlet state 
P= {*){% 

E x {a) = E 2 (P) =0, (17) 
E(a,f3) = -cos2(a-/3), (18) 

for which max^a'^^/ S(a, a', /3, /3') = 2^2, confirming 
that the singlet is a quantum state with maximal entan- 
glement. 

Analysis of the experimental data accor ding to the pro- 
cedure sketched earlier 0, HI, HI M Ep, Hi H|, yields 
results that are in good agreement with Ei(a) = E 2 (P) = 
and E(a,(3) = — cos2(a — (3), leading to the conclusion 
that in a quantum theoretical description, the density 
matrix does not factorize, in spite of the fact that the 
photons are spatially and temporally separated and do 
not interact. 



D. Classical simulation model 

A concrete simulation model of the EPRB experiment 
sketched in Fig. [¥] requires a specification of the infor- 
mation carried by the particles, of the algorithm that 
simulates the source and the observation stations, and 
of the procedure to analyze the data. In the following, 
we describe a slightly modified version of the algorithm 
proposed in Ref. [l3| . tailored to the case of photon po- 
larization. 

Source and particles: The source emits particles 
that carry a vector S ny i = (cos(£„ + (i — l)7r/2), sin(£ n + 



(i — l)7r/2), representing the polarization of the photons 
that travel to station i = 1 and station i = 2, respectively. 
Note that S n ,i 'S nj 2 = 0, indicating that the two particles 
have orthogonal polarizations. The "polarization state" 
of a particle is completely characterized by £„, which is 
distributed uniformly over the whole interval [0, 27r[. For 
this purpose, to mimic the apparent unpredictability of 
the experimental data, we use uniform random numbers. 
However, from the description of the algorithm, it will 
be clear that the use of random numbers is not essen- 
tial. Simple counters that sample the intervals [0, 2n[ in 
a systematic, but uniform, manner might be employed as 
well. 

Observation station: The electro-optic modula- 
tor in station i rotates S nj i by an angle jn^, that is 
a„ jt = (cos7 H) i, sinjn^). The number M of different 
rotation angles is chosen prior to the data collection 
(in the experiment of Weihs et ai, M — 2 [21]). We 
use 2M random numbers to fill the arrays (ai, cum) 
and (Pi, Pm)- During the measurement process we 
use two uniform random numbers 1 < m, ml < M 
to select the rotation angles 7„ i — a m and 7 n> 2 = 
(3 m i. The electro-optic modulator then rotates S n ,i — 
(cos(£„ + (i-l)7r/2),sin(£„ + (i-l)7r/2) by 7„ j4 , yielding 
S„,i = (cos(£„-7„ i4 +(i-l)7r/2),sin(£„-7„,i-|-(i-l)7r/2). 

The polarizer at station i projects the rotated vector 
onto its x-axis: S„^ • = cos(£„ — j n> i + (i — l)7r/2), 
where x.; denotes the unit vector along the x-axis of the 
polarizer. For the polarizing beam splitter, we consider 
a simple model: If cos 2 (£„ — 7„^ + (i — 1)tt/2) > 1/2 
the particle causes -D+i,i to fire, otherwise D-\^ fires. 
Thus, the detection of the particles generates the data 
x„,i = sign(cos2(£„ - 7^ + (i - 1)tt/2)). 

Time-tag model: To assign a time-tag to each event, 
we assume that as a particle passes through the detection 
system, it may experience a time delay. In our model, the 
time delay t n ^ for a particle is assumed to be distributed 
uniformly over the interval [to, io+T] , an assumption that 
is not in conflict with available data [27] . In practice, 
we use uniform random numbers to generate t„.i. As in 
the case of the angles £„, the random choice of t n> i is 
merely convenient, not essential. From Eq.©, it follows 
that only differences of time delays matter. Hence, we 
may put to = 0. The time-tag for the event n is then 
t n ,i £ [0,T]. 

There are not many options to make a reasonable 
choice for T. Assuming that the particle "knows" its 
own direction and that of the polarizer only, we can con- 
struct one number that depends on the relative angle: 
S n> i ■ Xi. Thus, T = T(£„ - 7^) depends on - j nii 
only. Furthermore, consistency with classical electro- 
dynamics requires that functions that depend on the 
polarization have period ir [TtJ- Thus, we must have 
T(£„ - 7r« + (i ~ 1)tt/2) = F{{S n>i ■ x,) 2 ). We already 
used cos2(£ n — 7^ + (i — l)7r/2) to determine whether 
the particle generates a +1 or —1 signal. By trial and er- 
ror, we found that T(£„ - 0i) = T F(|sin2{f„ — 6^)11 = 
T 1 sin 2 (£„ - 6»i ) | d yields useful results [H, [lj, GJ, M, Hi . 
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FIG. 8: (color online) Comparison between computer simula- 
tion data (red bullets) and quantum theory (black solid line) 
for the two-particle correlation E(a,/3). 



two different random number generators were used to se- 
lect the angles a m and (3 m >. The statistical correlation 
between m and m' was measured to be less than 10 -6 . 

From Fig. [51 it is clear that the simulation data for 
E(a, (3) are in excellent agreement with quantum the- 
ory. Within the statistical noise, the simulation data 
(not shown) for the single-spin expectation values also 
reproduce the results of quantum theory. 

Additional simulation results (not shown) demonstrate 
that the kind of models described earlier are capable of 
reproducing all the results of quantum theory for a sys- 
tem of two S=l/2 particles [jj, Q EE HI |37|. Fur- 
thermore, for W — > and in the limit that the number 
of events goes to infinity, one can prove rigorously that 
these simulation models give the same expressions for the 
single- and two-particle averages as those obtained from 
quantum theory 0, G3, Q3, II • 



F. Discussion 



Here, To = maxg T(9) is the maximum time delay and de- 
fines the unit of time, used in the simulation and d is a 
free parameter of the model. In our numerical work, we 
set T = 1. 

Data analysis: For fixed N and M, the algorithm 
generates the data sets just as experiment does [2l| . 
In order to count the coincidences, we choose a time-tag 
resolution < r < To and a coincidence window r <W . 
We set the correlation counts C xy (a m , (3 m >) to zero for 
all x, y — ±1 and m, m' = 1, M. We compute the dis- 
cretized time tags k n ^ = \t n ^/T~\ for all events in both 
data sets. Here \x] denotes the smallest integer that is 
larger or equal to x, that is \x~\ — 1 < x < \x~\ . Ac- 
cording to the procedure adopted in the experiment [2l| , 
an entangled photon pair is observed if and only if 
\k n ,i - k n ,2\ < k = \W/t\. Thus, if |fc„,i - fc„, 2 | < k, 
we increment the count C Xn x >Xn _ 2 (a m , /3 TO / ) . 

E. Simulation results 

The simulation proceeds in the same way as the exper- 
iment, that is we first collect the data sets {Ti, T2}, and 
then compute the coincidences Eq. ([SJ and the correla- 
tion Eq. (fTTj) . The simulation results for the coincidences 
C xy (a,(3) depend on the time-tag resolution r, the time 
window W and the number of events N, just as in real 
experiments [U M, H, 11 H HE M, H| . 

Figure [5] shows simulation data for E(a, (3) as obtained 
for N = 10 6 and W = r = 0.00025T . In the experiment, 
for each event, the random numbers A n ^ = 1, . . . , M se- 
lect one out of four pairs {(a*, (3j)\i, j = 1, M}, where the 
angles an and (3i are fixed before the data is recorded. 
The data shown has been obtained by allowing for M — 
20 different angles per station. Hence, forty random num- 
bers from the interval [0,360[ were used to fill the arrays 
(ax, ■ ■ ■ , um) and (j3i, . . . , (3m)- For each of the N events, 



Starting from the factual observation that experimen- 
tal realizations of the EPRB experiment produce the data 
{Ti,T2} (see Eq. ([7])) and that coincidence in time is a 
key ingredient for the data analysis, we have described a 
computer simulation model that satisfies Einstein's crite- 
rion of local causality and, exactly reproduces the correla- 
tion E(a.i, a.2) = —&i a.2 that is characteristic for a quan- 
tum system in the singlet state. Salient features of these 
models are that they generate the data set Eq. (J7J) event- 
by-event, use integer arithmetic and elementary mathe- 
matics to analyze the data, do not rely on concepts of 
probability and quantum theory, and provide a simple, 
rational and realistic picture of the mechanism that yields 
correlations such as Eq. (fT8)) . 

We have shown that whether or not these simulation 
models produce quantum correlations depends on the 
data analysis procedure that is performed (long) after 
the data has been collected: In order to observe the cor- 
relations of the singlet state, the resolution r of the de- 
vices that generate the time-tags and the time window 
W should be made as small as possible. Disregarding 
the time-tag data (d = or W > T ) yields results that 
disagree with quantum theory but agree with the models 
considered by Bell [26] . Our analysis of real experimen- 
tal data and our simulation results show that increasing 
the time window ch ang es the nature of the two-particle 
correlations [H Q GSM H3 ■ 

According to the folklore about Bell's theorem, a pro- 
cedure such as the one that we described should not ex- 
ist. Bell's theorem states that any local, hidden variable 
model will produce results that are in conflict with the 
quantum theory of a system of two S = 1/2 particles f26| . 
However, it is often overlooked that this statement can 
be proven for a (very) restricted class of probabilistic 
models only. Indeed, minor modifications to the orig- 
inal model of Bell lead to the conclusion that there is 
no conflict [H, El|. In fact, Bell's theorem does not 
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necessarily apply to the systems that we are interested 
in as both simulation algorithms and actual data do not 
need to satisfy the (hidden) conditions under which Bell's 
theorem hold [HSl, Si| . 

The apparent conflict between the fact that there exist 
event-based simulation models that satisfy Einstein's cri- 
terion of local causality and reproduce all the results of 
the quantum theory of a system of two 5=1/2 particles 
and the folklore about Bell's theorem, stating that such 
models are not supposed to exist dissolves immediately if 
one recognizes that Bell's extension of Einstein's concept 
of locality to the domain of probabilistic theories relies 
on the hidden, fundamental assumption that the absence 
of a causal influence implies logical independence [45j |. 
Indeed, in an attempt to extend Einstein's concept of a 
locally causal theory to probabilistic theories, Bell im- 
plicitly assumed that the absence of causal influence im- 
plies logical independence. In general, this assumption 
prohibits the consistent application of probability theory 
and leads to all kinds of logical paradoxes [H |47[ ■ How- 
ever, if we limit our thinking to the domain of quantum 
physics, the violation of the Bell inequalities by experi- 
mental data should be taken as a strong signal that it is 
the correctness of this assumption that one should ques- 
tion. Thus, we are left with two options: 

• One accepts the assumption that the absence of a 
causal influence implies logical independence and 
lives with the logical paradoxes that this assump- 
tion creates. 

• One recognizes that logical independence and the 
absence of a causal influence are different con- 
cepts [H |4j 13 and one searches for rational ex- 



planations of experimental facts that are logically 
consistent, as we did in our simulational approach. 

IV. CONCLUSION 

The simulation models that I described in this talk 
do not rely on concepts of probability theory or quan- 
tum theory: They are purely ontological models of quan- 
tum phenomena. The salient features of these simulation 
models @, B [1 El El El [M IH are that they 

1. generate, event- by-event, the same type of data as 
recorded in experiment, 

2. analyze data according to the procedure used in 
experiment, 

3. satisfy Einstein's criterion of local causality, 

4. do not rely on any concept of quantum theory or 
probability theory, 

5. reproduce the averages that we compute from 
quantum theory, 

We may therefore conclude that this computational mod- 
eling approach opens new routes to ontological descrip- 
tions of microscopic phenomena. 
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